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Abstract 

In this letter we extent the overrelaxation algorithm, known to be very efficient in 
gauge theories, to coupled gauge-Higgs systems with a particular emphasis on the update 
of the radial mode of the Higgs field. Our numerical tests of the algorithm show that the 
autocorrelation times can be reduced substantially. 
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The electroweak baryogenesis has attracted a lot of attention recently [l|] . In order to understand 
this phenomenon a detailed description of the electroweak phase transition (PT) is needed. In 
its study it is sufficient to concentrate on the SU(2) group alone and to neglect all fermion 
fields and the U(l) subgroup. One therefore remains with a coupled system of gauge and 
Higgs fields which are elements of SU(2) and complex doublets, respectively. Although this 
so-called SU(2)-Higgs model has been investigated up to 2-loop order at finite temperature in 
perturbation theory by now ||, numerical simulations, which intrinsically contain also non- 
perturbative effects, are desirable. This is even more so as the 2-loop results showed large 
corrections to quantities like the surface tension as compared to 1-loop results. Therefore it 
seems necessary to confront the perturbative approach with numerical "experiments" to test 
the reliability of the perturbative results. 

The Monte-Carlo simulations are necessarily done on a discretized (euclidean) space time 
lattice. To obtain results which are not distorted by the lattice one would like to go close to 
the continuum and work in the "scaling region". To approach this limit some of the correlation 
lengths £j have to diverge. This causes the basic problem of all numerical simulations as with 
growing correlation lengths the problem of critical slowing down arises. In order to obtain 
independent configurations, suitable for measurements of physical quantities, one has to take 
into account the autocorrelation time which is the number of iterations with a given algorithm 
to reach a new independent configuration. The autocorrelation time grows with the correlation 
length as r oc f with z the so called dynamical critical exponent. For local algorithms like 
Metropolis or heatbath z is known to be 2. This implies that for £'s of the order of 5-10, 
which are realistic values of todays simulations of the SU(2)-Higgs model, the autocorrelation 
time can be of the order of 25-100. Indeed, as will be shown below, in the SU(2)-Higgs model, 
for parameter values where one can compare results from perturbation theory and numerical 
simulations, the autocorrelation times are 0(100). This appears to be a major drawback for 
Monte-Carlo studies of the electroweak PT. 

It is consequently not surprising that fighting critical slowing down is one of the major 
activities in the area of Lattice field theory. In fact, in the last few years several important 
steps have been done to solve this problem. Cluster |J| and multigrid M techniques are able 
to reduce critical slowing down almost completely, giving z ~ 0. However, these techniques are 
so far either only applicable for spin models (cluster) or did not lead to a big improvement in 
non-abelian gauge theories (multigrid). 

For gauge theories another interesting approach, the overrelaxation method, as initiated 
by Adler || as a generalization of the heatbath algorithm seems to be most promising. In its 
popular limiting case the field evolution is deterministic and conserves energy (microcanonical). 



In order to restore the necessary ergodicity one uses standard, ergodic updates (e.g. Metropolis 
or heatbath) and microcanonical updates alternately with a given mixing ratio. These so called 
hybrid algorithms turned out to be extremely efficient in pure gauge theories 0. In this case 
the overrelaxation step is a kind of reflection of the local field with the following features: 

a) the energy is unchanged 

b) since the reflection is an element of the gauge group, the measure in the defining functional 
integral is invariant 



We will extent this approach to coupled gauge-Higgs systems. Consider an SU(2) gauge-Higgs 
system with the following lattice action: 

S = Sg + Sh, (1) 

where 

S B = -~T±(ZU P ) (2) 
z v 

is the usual Wilson plaquette action and 

S h = -KTr(]T $tU Xj ^ x+ ») + £ [*t$ a + A ($t $ x - I) 2 ] (3) 

X,fl X 

describes the self interaction of the four component Higgs field $ x and its coupling to the gauge 
field. Here one can introduce the standard 2x2 matrix notation, $> x = p x a x , where p x e R, 
p x > 0, and a x G SU{2). 

The reflection of the gauge field and the angular part of the Higgs field can be done anal- 
ogously as for pure gauge fields and the above a) and b) conditions can be satisfied. However, 
for the radial part of the scalar field the reflection fulfills a) but not b), thus a careful treatment 
of the measure is needed in order to obtain a proper technique. 

First we will study the microcanonical updating of the link variables. The action can be 
written as a sum of two terms: a term which contains the local link variable and a constant 
term. 

S = TT{U Xtll V) + const, (4) 

where 

V = ~2 U x+^,uU ] x+u ^Ul u - kJ2 Pxpx+^a x+ll al. (5) 

Consider an update as a "reflection" of U 

= vM*vl (6) 
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where V is the normalized V, thus qV = V, where q G R and V G SU(2). (Note the special 
feature of the SU(2) group that the sum of several SU (2) matrices with real coefficients is a 
product of a real number and an SU(2) matrix.) As in the case of pure gauge theory it is easy 
to see that this updating satisfies the above mentioned a) and b) conditions. 

The micro canonical updating of the angular part of the scalar field can be done completely 
analogously. 

S = Tr(a x V) + const, (7) 

where 

V = -np x Px+^U x ^a x+fM , (8) 

thus the updated angle 

ol x = VoalVo, (9) 

where Vq is, as above, the normalized V . Again, one can show that the above a) and b) 
conditions are satisfied. 

The most interesting case is the microcanonical updating of the radial mode of the scalar 
field. Separating the angular and the radial modes as $ = pa, the measure will get a 

form of dadpp 3 . Therefore one gets a potential 

V{ Px ) = -Cp x + p 2 x + \{pl - l) 2 - 3log(p x ), (10) 

where 

c = /^E^Xm^) • ( n ) 

We have plotted this local V(p) potential in fig. I. for A = 0.0001 and for a typical 
C = 2 value. The potential has one minimum, but is not symmetric with respect to its 
minimum. By a simple analysis of V(p) one can show, that it has only one minimum for 
A < (3/2 + y/Tl/2) 2 ps 9.98. For even larger values of A, one needs k ps (9(10) to have a second 
minimum. These k-s are clearly out of any range of physical interest. Thus, we can assume 
that V(p) has only one minimum. 

The procedure to perform the overrelaxation step in p is now clear. For each p x calculate the 
value of p' x 7^ p x such that V(p x ) = V(p' x ). Note that this is not just a reflection with respect 
to the minimum of the potential as V(p) is not symmetric. In practice we found that a reliable 
way to determine p is to calculate the first, second and third derivatives of V (analytically) 
and calculate p' x for which the Taylor polynom gives the same value for the potential. After 
this starting step one-two Newton iterations give the proper p' x value with a relative accuracy 
of 10~ 6 . If the first step can not be made, or turns out to be too rude Newton alone or bisection 
can be used. 

The non-trivial part of the overrelaxation step in p is the observation that p is a nonlinear 
function of p. Therefore, the transformation of p to p does not leave the measure invariant. 
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As a consequence, one has to correct for this by a reject/accept step. In order to obtain the 
proper updating procedure let us consider the stability equation for p. 



JdpT(p,p)P(p)=P(p), (12) 

where P(p) is the equilibrium distribution in p and T(p, p) is the overrelaxation operator. It 
has the form 

where M{p) is the "mirror image" of p, thus 

V(M(p)) = V(p) , but p^M(p), (14) 

unless p is exactly at the minimum, where p = M(p), of course. Evaluating the above stability 
integral on gets 

\ipMMp - m(p))p(p) = m^^m- (15) 

Combining this with the stability equation one obtains 

A(, ) = |tf W | = |->| = |_ / _|. (16) 

The Monte-Carlo realization of the above T(p,p) overrelaxation operator can be done in the 
following way. One determines M(p) "mirror image" of p and calculates the V' (p) and V' (M(p)) 
derivatives. The p — > M(p) updating is accepted if 

,17) 



•V(M(p)) 



where r is a random number between and 1. Otherwise the updating is rejected. In practice 
we found the acceptance always to be larger than 80%. 

Note, that the same result can be obtained by considering the probability densities. The 
probability of finding the system in a state with local energy between V and V + AV is a sum 
of two probabilities. The probability to find it between p and p + Ap plus the probability to 
find it between M(p) and M(p) + AM(p), where V{p) = V{M(p)) = V and V(p + Ap) = 
V(M(p) + AM(p)) — V + AV. Since we want the algorithm to be microcanonical, dV = 
V' \p)dp — V' (p')dp' — (p — M{p)). Clearly one regains (16) for A(p) and the reject/accept 
step (17). 
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We have tested the algorithm described above in the SU(2)-Higgs model. For parameter values 
A = 0.0001, (3 = 8 and k = 0.129 the system is in the Higgs region of the model. We measured 
the normalized autocorrelation function Tq of some operator O which is defined as 

< O(0)O(t) > - < O(0) > 2 

To{t) < o(o) >2 (18) 

where t indicates a fictitious Monte Carlo time corresponding to the number of sweeps. Tq falls 
off exponentially 

r (t) oc exp(-t/r) for t -» oo (19) 

which defines the exponential autocorrelation time r. In fig. 2 we plot the autocorrelation 
function T p for the length of the Higgs field p. This operator in general shows the longest 
autocorrelation time in the SU(2)-Higgs model. The solid line corresponds to a pure Metropolis 
simulation and the dashed line shows the result of a hybrid overrelaxation with a mixing ratio 
of Metropolis to overrelaxation 1:1. The measurements have been performed after each sweep 
of a Metropolis or an overrelaxation update through the lattice. Although the lattice used is 
quite small (8 4 ) the autocorrelation function for the Metropolis algorithm shows a very slow 
fall off, indicating an autocorrelation time of about TMet ~ 130. On the other hand the hybrid 
overrelaxation algorithm gives a considerable improvement, tor ~ 25, over the Metropolis 
algorithm |. The results for the exponential autocorrelation times r eq.(19), given above, have 
been checked against the integrated autocorrelation times and a complete agreement has been 
found. 

The dashed curve in fig. 2 corresponds to a combination of only one Metropolis to one over- 
relaxation step. We want to note that the main effect of the improvement stems from the 
overrelaxation in p alone. Switching off the overrelaxation in U and a we found an almost 
identical autocorrelation function. This clearly demonstrates that the operator p gives the 
slowest mode. Increasing the mixing ratio of Metropolis to overrelaxation steps does not im- 
prove the autocorrelation time. This is understandable from the p-potential eq.(10) as more 
overrelaxation steps reflects p only back and forth, not leading to substantial changes in the 
configuration. 



In simulations of the SU(2)-Higgs model, which is the most relevant part for studies of the 
electroweak phase transition, using a Metropolis algorithm one detects large autocorrelation 



2 Based on the ideas of the present letter, a hybrid of heatbath and overrelaxation was analyzed in || and a 
similar improvement was found. 



times already on small lattices. In order to be able to have reliable numerical results it is 
therefore necessary to find an algorithm which improves this behaviour. We have extended 
the overrelaxation algorithm which seems to work very efficient in pure gauge theories to the 
SU(2)-Higgs model. The treatment of the radial part of the Higgs-field in the algorithm led to 
an additional accept/reject step which corrects for the measure term. 

We have tested our version of the overrelaxation algorithm in the SU(2)-Higgs model at 
parameter values interesting for comparisons with perturbation theory. As fig. 2 shows, the 
autocorrelation time can be reduced substantially, namely from r ~ 130 in the Metropolis 
algorithm (solid line) to r « 25 in the overrelaxation algorithm. We want to emphasize that 
this result is obtained on a small (8 4 ) lattice. It is expected that the hybrid overrelaxation 
algorithm has a dynamical critical exponent of z ~ 1 as compared to the Metropolis or heatbath 
with z«2. Therefore one can expect even better improvements on larger lattices. 

We found it sufficient to perform only one overrelaxation step in the length of the Higgs field 
p to obtain the best improvement. We expect, however, that as in pure gauge theories || |3|, 
one can obtain even smaller autocorrelation times by tuning the mixing ratio of Metropolis to 
overrelaxation steps for the update of the gauge field U and the angle of the Higgs field a. This 
tuning will depend on the parameter values where simulations are performed. The optimal 
mixing ratio should be determined from case to case. 

We think that our numerical results are very promising and that the hybrid overrelaxation 
we are suggesting is a big improvement over algorithms used so far for simulations of the SU(2)- 
Higgs model. We hope that by careful fine tuning of the mixing ratio of different algorithms 
for the update of the variables U, p and a the autocorrelation times might be reduced to small 
numbers r < 10 as it can be obtained in pure gauge theories and one can therefore obtain 
high precision numerical data to shed new light on the the electroweak phase transition beyond 
perturbation theory. 
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Figure Caption 

Fig.l The local potential eq.(10) as a function of the radial length of the scalar field p for 
A = 0.0001 and C = 2. 

Fig. 2 The logarithm of the autocorrelation function for the length of the Higgs field p as a 
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function of the Monte Carlo "time" t. The parameter values are (3 — 8, A = 0.0001 and 
k = 0.129 and belong to a point in the Higgs region of the model. The lattice is 8 4 . The solid 
line is the Metropolis algorithm alone. The dashed line is a hybrid of one Metropolis to one 
overrelaxation step. 
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